Mon. Not. R. Astron. Soc. 000, 1-12 (0000) Printed 2 February 2008 (MN WF$L style file vl.4) 



Active galactic nuclei and the minor merger hypothesis 



Philip Kendall 1 , John Magorrian 1 ' 2 and J.E. Pringle 

1 Institute of Astronomy, Madingley Road, Cambridge, CB3 OH A, UK 

2 Department of Physics, University of Durham, South Road, Durham DH1 3LE 



Accepted — . Received — ; in original form — . 



ABSTRACT 

We have investigated the dynamics of the merging process in the minor merger 
hypothesis for active galactic nuclei. We find that for a satellite galaxy to be able to 
merge directly with the nucleus of the host galaxy (for example, to give rise to the 
compact dust discs which are seen in early type active galaxies) requires the initial 
orbit of the satellite to be well aimed. For the case of the host galaxy being a disc 
galaxy, if the initial orbits of the satellites are randomly oriented with respect to 
the host galaxy then the orbits of those which reach the host nuclear regions in a 
reasonable time, are also fairly randomly oriented once they reach the nucleus. We 
note that this result might be able to provide an explanation of why the jet directions 
in the nuclei of Seyfert galaxies are apparently unrelated to the plane of the galaxy 
discs. 

Key words: galaxies: interactions, galaxies: individual: NGC 3379, hydrodynamics, 
methods: numerical 



1 INTRODUCTION 

In this paper we investigate the hypothesis that activity 
in galactic nuclei is induced by the acquisition and dis- 
ruption of a small gas-rich companion galaxy. There has 
been considerable discussion of this 'minor merger' hypoth- 
■ esis in the literature (see, for example, recent discussions 
in Taniguchi, 1999; Chatzichristou, 2000, a,b; 2001 a,b, and 
papers therein). 

In elliptical galaxies, a study of the effects of minor 
mergers is motivated by the observational evidence (van 
Dokkum & Franx, 1995) that dust discs are present in the 
cores of a high fraction of ellipticals, with a higher detec- 
tion rate in radio galaxies (Verdoes Kleijn et al., 2000). This 
result is consistent with the results of surveys of 3CR ra- 
dio galaxies by de Koff et al. (2000) and by Martel et al. 
(1999). The simplest interpretation of these is that the dust 
structure seen in these galaxies is the debris trail of a small 
merger incident. The analysis of de Koff et al. (2000), and 
the finding by Sparks et al. (2000) that galaxies with opti- 
cal jets tend to display face-on (i.e. round) dust discs, both 
suggest that the radio jets in radio galaxies show a tendency 
to be aligned perpendicular to the observed dust disc struc- 
ture. In this picture a recent merger would have stimulated 
the current nuclear activity. The surrounding dust structure 
would indicate the plane of the orbit of the disrupted satel- 
lite, and the resulting radio jets would naturally be expected 
to be perpendicular to this plane. 



In view of this, we investigate first the conditions re- 
quired for a small gas-rich satellite galaxy to be accreted by 
a larger galaxy in such a way that it gives rise the a dust disc 
similar to those observed. To focus the discussion we con- 
centrate on modelling the formation of the dust disc in the 
otherwise unremarkable El galaxy NGC 3379. Our results 
are straightforwardly applicable to other gas-poor ellipticals. 
The dust disc in NGC 3379 (van Dokkum & Franx, 1995; 
Ferrari et al., 1999) has a radius ~ 100 pc, and a mass esti- 
mated to be around 150 Mq. In addition Ferrari et al. (1999) 
show that there is more patchy dust extending out to about 
1 kpc from the nucleus. We adopt a simplified procedure to 
model the interaction. We use spherical models for both the 
major galaxy and the minor mergee and model the motion 
of the orbit of the satellite using simple dynamical friction, 
along with a straightforward model for the tidal disruption 
of the satellite (Section 2). We use a simple model for the 
subsequent evolution of the stripped gas (Section 3). Our 
basic conclusion (Sections 5.1 and 6) is that in order to gen- 
erate a dust disc of the kind observed it is necessary for the 
initial trajectory of the merging satellite to be accurately 
directed towards the nucleus of the major galaxy. 

However, reality turns out to be more complicated than 
this simple picture might suggest. First there is the recent 
analysis by Schmitt et al. (2002) who consider only those 
E/S0 galaxies with active nuclei which have well — defined 
dust discs in their central regions, rather than galaxies whose 
dust distribution is more dispersed. Schmitt et al. conclude 
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that the radio jets from the centres of these discs are not in 
fact oriented perpendicular to the dust discs' axes, and the 
angles the jets make with the axes are consistent with being 
distributed evenly in the range 0° — 75° . Second there is the 
finding by Kinney et al (2000) that in Seyfert galaxies the di- 
rections of the central jets are consistent with being oriented 
randomly in space, independent of the plane of the gas disc 
of the host spiral galaxy. If the minor merger hypothesis is 
to hold, then the apparent misalignments between between 
the angular momentum of the accreting material (as evi- 
denced by the observed discs in early type galaxies and by 
the disc plane in Seyferts) and the observed direction of the 
jets needs some explanation. One possible explanation is the 
idea that the jet direction is governed by the spin of the cen- 
tral black hole (through the Bardeen-Petterson effect (e.g. 
Scheuer and Feiler, 1996 and references therein) and that 
for some reason the spin of the black hole is misaligned with 
the current accretion flow. While this is a reasonable expla- 
nation if the galaxy is spherically symmetric it would seem 
difficult to achieve if the spin of the black hole has been gen- 
erated by the gas accretion process and if the potential of 
the galaxy is such the the accreted gas would (in the long 
term) be expected to have a preferred direction for its angu- 
lar momentum vector. One way round this problem, follow- 
ing the ideas of Wilson & Colbert (1995; see also Merritt, 
2002), is if the black hole spin is generated not by accretion 
of gas, but rather by accretion of black holes presumably 
contained in the nuclei of the small merging satellites. Even 
a small black hole (say a tenth of the mass of that in the 
host galaxy) can in its final approach add substantial orbital 
angular momentum to the spin of the final coalesced black 
hole. For this process to work, what is required is that the 
central (nuclear) black hole of a merging satellite be able to 
reach the nuclear regions of the major galaxy with the an- 
gle between its orbital angular momentum and the axis of 
symmetry of the gravitational potential of the major galaxy 
distributed over a wide range of angles. 



To investigate this possibility we extend our analysis 
to the case where the major host galaxy is a spiral galaxy. 
Here we are not concerned with the fate of any gas that 
might have been in the merging satellite, but are, rather, 
concerned with the fate of its nuclear black hole. We make 
use of similar simplified dynamics using the concepts of grav- 
itational drag and a similar approximation to take account 
of tidal stripping to those we used before, but we now ap- 
ply these concepts to a a galaxy with a more complicated 
internal structure (Section 4). As before we find that the 
cross-section for an incoming satellite to reach the nucleus 
of the host galaxy in a Hubble time is relatively low. How- 
ever here we pose the additional question: of those satellite 
nuclei which do reach the central regions of the host galaxy 
in a reasonable time, what is the distribution of the incli- 
nations of their orbits relative to normal to the disc of the 
galaxy? In order for the merger hypothesis to remain viable 
we shall require the distribution we find to very broad. We 
summarise our findings and give discussion of our results in 
Sections 5.2 and 6 respectively. 



2 MERGING WITH A SPHERICAL GALAXY 

In this Section we describe our simplified galaxy dynamics 
and merging procedure. 



2.1 Galactic models 

As mentioned in Section 1, we first investigate the particular 
case in which the major galaxy is spherically symmetric as a 
reasonable approximate model for the merging of a satellite 
with the galaxy NGC 3379. 



2.1.1 Major galaxy (NGC 3379) 

To model the density profile of NGC 3379 we have made 
use of Lima Neto, Gerbal and Marquez's (1999) analytic 
approximation to the deprojection of its surface brightness 
profile obtained with HST (Pastoriza et al. 2000)*. Lima 
Neto et al. (1999) assume that the galaxy is spherical, with 
a constant mass-to-light ratio T, and find that 



p(r)=p (£) "expf-QJ 



(1) 



where v — 0.424 and p — 0.572. Using the normalization 
to the surface brightness given by Pastoriza et al. (2000), 
together with an assumed distance to NGC 3379 of 9.9 Mpc 
and a mass to light ratio of T = 5.3T0 (Magorrian et al. 
1998), gives Po = 3.96 x lO 3 M /pc 3 , and a = 18.6 pc. The 
resulting visible mass is 2.9 x 1O 1O M , with a half- mass 
radius of about 800 pc, excluding any dark halo. We find that 
the inclusion of a dark halo or a central 10 s Mq black hole 
(Magorrian et al. 1998; Gebhardt et al. 2000) has negligible 
effect on our results below. 

We calculate the mass distribution, potential and veloc- 
ity dispersion of the galaxy through the usual equations 
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where we have assumed for simplicity that the galaxy's ve- 
locity dispersion tensor is isotropic. We evaluate the inte- 
grals adaptively using Simpson's rule until an accuracy of 
one part in 10 6 is reached and store the calculated values 
of M(r), $(r) and a(r) on grids with 30 points spaced log- 
arithmically in radius between 100 pc and 100 kpc. Values 
between tabulated points are obtained via cubic spline in- 
terpolation. 



2.1.2 Satellite 

We have considered two different models for the small in- 
falling satellite galaxy. One is a rigid Plummer sphere with 
density distribution given by 



* Note that there is a typographical error in Pastoriza et al. 
(2000): in equation 2 'cxp' should be '10 to the power of (Ferrari, 
private communication) 
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and the other is a rigid Hernquist sphere with density dis 
tribution given by 

M a 1 



p(r) = 



2tt r ( a + r) 



(6) 



For both of these cases we calculate the mass and velocity 
dispersion as for the host galaxy, except with a grid spaced 
logarithmically between 0.1 kpc and lOkpc. The total mass 
(M = 10 9 M Q ) and scale lengths (b = 0.4 kpc and a = 
0.23 kpc) were chosen to match the estimated mass and half- 
mass radius (0.55 kpc) of Sagittarius (Helmi & White 2000), 
as a typical dwarf galaxy. We will see in Section 2.3 that the 
choice of model for the satellite does not make a significant 
difference to our results. 



2.2 The orbit of the infalling satellite 

In computing the orbit of the infalling satellite, we neglect 
changes which might be induced by the satellite in the host 
galaxy. This approximation has been the subject of debate 
in the literature (Zaritsky & White 1988; Hernquist & Wein- 
berg 1989), with the conclusion being reached that it is rea- 
sonable to ignore the galaxy response for low mass satellites 
(Velazquez & White 1999). Physically, this is what is ex- 
pected, as the response induced in the galaxy is proportional 
to the mass ratio of the satellite to the galaxy M 3 /M g , and 
the effect of this response is again proportional to the mass 
ratio, so the overall effect will be proportional to (M s /M g ) 2 . 

We model the motion of the satellite using Chan- 
drasekhar's dynamical friction formula (Binney & Tremaine 
1987). Thus 

dv B _ kM s p 



(7) 



where v s is the velocity of the satellite, p is the density of 
the medium, k is given by 



(8) 



k = 47T In AG 2 (&i{X) - ^| exp (-X 2 )^j , 

G is Newton's gravitational constant, X, the dimensionless 
speed of the satellite, is given by 



X = 



(9) 



and A is the Coulomb logarithm discussed below. 

Equation (7) applies to the case of a point mass 
moving through an infinite homogeneous medium with a 
Maxwellian velocity distribution. However, numerical exper- 
iments (Velazquez & White 1999) show that it also provides 
a reasonable approximation to the gravitational drag expe- 
rienced in more realistic situations, provided one chooses the 
Columb logarithm In A appropriately. To determine whether 
it is adequate for the situation we are considering, we use 
it to reproduce the results from the infall of a satellite into 
a dark halo when the galaxy was modelled using a full N- 
body treatment (van den Bosch et al. 1999). We note that 
this procedure has been investigated in some detail by Tay- 
lor & Babul (2001), who demonstrate that it does provide a 
reasonable approximation. 
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Figure 1. A comparison our approximate method for follow- 
ing the effects of dynamical friction on a satellite's orbit (right) 
against results from a full TV-body simulation (van den Bosch et 
al. 1999; left panels). The top panels show how the distance of 
the satellite from the centre of the host varies with time. The 
bottom two show the eccentricity of the satellite's orbit, defined 
as e = (r+ — r_)/(r+ +r_), where r+ and r_ are the satellite's 
apo- and peri-centre radii respectively. 



To perform the integration, we used the algorithm from 
the fifth-order Cash-Karp Runge-Kutta integrator 'odeint' 
(Press et al. 1992). The timestep used is automatically ad- 
justed by the routine to control errors, and in no case did 
the timestep reach a significant fraction of the dynamical 
time of the satellite. 

As a test of our procedure, we use our method to model 
van den Bosch et al.'s (1999) Plummer sphere sinking into 
a truncated isothermal halo. We find that we are able to 
reproduce their results over a wide range of conditions if we 
choose In A « 2. In Figure 1 we show a comparison between 
the full A-body calculations of van den Bosch et al (1999) 
and our approximate procedure for a typical case. As can 
be seen, there is reasonable agreement between the two. Us- 
ing these results we shall adopt the same choice In A = 2 for 
the orbital computations contained in this paper. We should 
note that, for the spherical case, our treatment of dynami- 
cal friction also ignores the effects of possible galaxy rotation 
and velocity anisotropy, but unless these are extreme we ex- 
pect that this approximation will be adequate in the current 
case. 



2.3 Tidal stripping 

In the model used by van den Bosch et al. (1999), the satel- 
lite maintained a constant mass all the way through the sim- 
ulation. Due to tidal effects, the outer regions of the satellite 
will not remain bound to the satellite. We are interested, 
therefore, in the manner in which the infalling satellite is 
tidally stripped as it passes through the host galaxy. This 
has an effect upon the orbital dynamics, because the decel- 
eration due to Chandrasekhar friction is proportional to the 
mass of the satellite. 

We use a simple algorithm to simulate the tidal strip- 
ping of the satellite. At each instant, we define the tidal 
radius, r t , of the satellite as being the radius at which the 
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Figure 2. Remaining satellite mass as a function of distance from 
the centre of the host galaxy obtained using the tidal stripping 
prescription (10). The solid and dashed curves plot the results 
for the Plummcr and Hcrnquist models for the satellite respec- 
tively (section 2.1.2). The former cuts off at around lkpc as its 
finite central density means that it is completely disrupted by 
that point. 



mean density of the satellite is equal to the mean density of 
the galaxy, viz., 



M s (r t ) _ M g (r 3 ) 



4 

— 7YT" 

3 



4 3 



(10) 



where M s (r) is the mass of the satellite contained within 
a radius r, M g (r) is the same for the galaxy and r s is the 
galactocentric radius of the centre of the satellite. We then 
truncate the satellite at r = r t , keeping the density structure 
internal to this point unchanged. This reduces the mass of 
the satellite to M s (r t ). In reality the satellite will react to 
this density truncation on an internal dynamical timescale, 
but since by construction^ this timescale is similar to the 
orbital timescale at this point, neglect of this effect is likely 
to be a reasonable approximation for current purposes. 

In Figure 2, we show the outcome of this stripping algo- 
rithm by plotting the fraction of the satellite galaxy which 
remains bound as a function of its distance from the centre 
of the host galaxy. We show the results for both models of 
the satellite galaxy. The Plummer sphere is completely dis- 
rupted by the time it comes within 1 kpc of the centre of 
the host galaxy. The Hernquist sphere loses 99% of its initial 
mass by the time it comes within is 400 pc from the centre 
of the galaxy. These distances are much greater than the 
size of the observed dust disc, which has a radius of around 
100 pc, but are not that dissimilar to the overall size of the 
dust spread within the galaxy. We note that these result de- 
pends only upon the density profiles of the satellite and the 



t The tidal radius is defined by equating the mean density of the 
satellite to the mean density of the main galaxy, and dynamical 
timescalcs just depend on mean densities. 



galaxy. They are just the result of solving Equation (10), 
and do not involve any consideration of the dynamics. 



3 DEPOSITION AND LOCATION OF THE 
STRIPPED GAS 

3.1 Low energy orbits 

To make a preliminary investigation of the initial conditions 
required for a merging satellite to be able to deposit signifi- 
cant gas into the nucleus of the host, we start by considering 
a series of low-energy orbits for the satellite galaxy. In each 
one the satellite starts on an orbit at an apocentre of 10 kpc 
from the galactic centre (that is, it starts with zero radial ve- 
locity at that point), and we vary the initial varying angular 
momentum. We characterize the initial angular momentum 
in terms of an eccentricity e, where e = (r+ — r_)/(r+ +r_) 
and r+(r_) are the apo- (peri-) centres of the orbit of a test 
particle with the same initial location and velocities. We 
use the models of NGC 3779 and the satellite as described 
in Section 2. 

To simplify the discussion we shall assume that the gas 
in the satellite galaxy is distributed spatially in the same 
way as the stellar matter. We make this assumption for two 
reasons: First, such an assumption probably overestimates 
the concentration of the gas within the dwarf. This will lead 
to the gas being stripped stripped from the dwarf later, and 
thus at a lower angular momentum than would otherwise 
occur. Hence this assumption will be an upper bound on the 
magnitude of any effect seen. As we show later, even with 
this assumption, we cannot replicate the observed effects. 
Second, by assuming a spherical dwarf, we can dramatically 
reduce the computational power required for our simulations 
and thus investigate a larger region of orbital parameter 
space. 

As an initial estimate of where the stripped gas ends 
up, the simplest assumption is that stripped material con- 
serves its angular momentum about the centre of the host 
galaxy (since the host is spherically symmetric), and that 
each element of gas, by losing energy through dissipative 
processes ends up in a circular orbit corresponding to its 
initial angular momentum L. The radius r c of this orbit is 
given by 



L 2 = GM g (r c )r c 



(11) 



The results of making this simple approximation are shown 
in Figure 3. The outermost spike in each of these mass dis- 
tributions is caused by the stripping which occurs at the 
first pericentre passage. For a high angular momentum or- 
bit (e = 0.5), the stripping occurs over a range of radii as 
the satellite's orbit gradually decreases in size due to gravi- 
tational drag. For the intermediate angular momentum orbit 
(e = 0.75) the stripping occurs more rapidly and thus the 
range of angular momentum in the stripped gas is reduced. 

In order to acheive the required radial penetration, we 
also investigated a low angular momentum orbit (e = 0.98) 
for which the satellite's orbit is almost radial. In this case, 
the gravitational drag plays almost no role, and the satellite 
is totally disrupted on the first periastron passage, but the 
gas still ends up at radii exceeding 250 pc. Thus we conclude 
that, using these approximations, the gas cannot end up 
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Figure 3. The surface density distribution of stripped gas if the gas is assumed to maintain its angular momentum, but loses energy 
until its orbit becomes circular. The panels show results for initial eccentricities e = 0.5 (left) and 0.75 (right). The large spike in every 
plot is due to the large amount of stripping which occurs during the first pericentre passage, whilst the spiky structure seen in the ring 
profiles is due to stripping occurring near successive pericentre passages. 



within around 100 pc unless the satellite was initially on a 
very eccentric, indeed essentially radial, orbit. 

3.2 Gas dynamics 

The assumption we made above, that an element of stripped 
gas, as it falls in, loses energy by interactions with other 
elements, but does not exchange angular momentum with 
them is not fully realistic. Since the angular momentum is 
crucial in determining the ultimate location of the gas, and 
since the discs we are obtaining are still much larger than 
observed discs we now consider a more realistic model for 
the radial evolution of the stripped gas. 

In the spirit of the approximations made in this paper, 
rather than undertaking a full hydrodynamical treatment 
(for example using SPH) , we opt for computational simplic- 
ity and have followed the approach of Lin and Pringle (1976) . 
In such a scheme, the gas particles move freely in the com- 
bined potential of the host galaxy and satellite, except that a 
form of dissipation is introduced which simulates thermally 
radiative shocks, by preventing the crossing of particle or- 
bits while exactly conserving mass and angular momentum. 
In this scheme, if two particles are close to each other, they 
interact, losing a fraction of their energy, but retaining their 
total momentum. One effect of this is to redistribute the 
angular momentum between different particles. For exam- 
ple, a circular ring of particles acts as an accretion disc and 
gradually spreads out to form a wider ring. In this paper, 
we consider inelastic collisions between particles in order to 
produce the greatest dissipation. This is clearly the most 
extreme situation, and therefore the results of this section 
will provide an upper bound on the magnitude of this effect. 

The scheme used to determine when collisions occur is 
simple: we divide the space in the host galaxy into cells, and, 
after every time step of the integration, if two (or more) par- 
ticles occupy the same cell, they interact. We use a cylindri- 
cal polar grid (7?, <f>, z) to describe the cells, with the z = 
plane aligned with the orbit of the satellite galaxy. The cell 
boundaries are spaced logarithmically in R, from 100 pc to 
20kpc with R n+ i/R n — 1.02, and linearly in the z from 
z — — 5.05kpc to z — 5.05 kpc with z n+ i — z n — 0.1 kpc. 
The choice of R n +i/R n is governed by a balance between 



the need for computational efficiency (if the zones are too 
small then there are two many cells with too few particles in 
them which permits interpenetration of fluid flows), and a 
need to keep the zones small, so that viscous effects do not 
become too rapid, resulting in unphysically rapid spread- 
ing of the gas discs which form. We performed runs with 
Rn+i/Rn of 1.05, 1.02 and 1.01. The results from the runs 
at 1.02 and 1.01 were very similar, implying that conver- 
gence had been achieved. Hence we choose R n +i/R n = 1-02 
to reduce the computational time required for our simula- 
tions. We choose 314 equispaced points in <j>. This means 
that the projection of the cells onto the z — plane is ap- 
proximately square and that the effective cross section of 
our gas particles varies from 2 pc at R — 100 pc to 400 pc at 
R = 20 kpc. 

We create one gas particle, of mass 10 5 Mq for every 
10° M Q stripped from the satellite. The particle is started 
with a velocity 

v g = v s +a s (\r g - r s \)f, (12) 

where v g and v 3 are the velocity of the gas particle and the 
satellite respectively, r g and r s are the positions of the cen- 
tres of the gas particle and the satellite, f is a unit vector 
in a random direction and a s (r) is the velocity dispersion 
of the satellite at a radius r. A particle is removed from the 
computation when one of the following conditions is met: 
the particle is within 100 pc of the centre of the galaxy (at 
which point, the gas is considered to have become part of the 
nuclear gas disc); the particle is within 90% of the current 
tidal radius of the satellite (when it is considered to have 
been reabsorbed by the satellite); or, the particle moves be- 
yond the outermost edges of the grid at R = 20 kpc and 
\z\ = 5.05 kpc. 

The radii at which the gas settles initially using this 
model are shown in Fig. 4. For comparison with the pre- 
vious estimates (Figure 3) we show the result for initial 
eccentricities of 0.5 and 0.75. We do not show the results 
for an initial eccentricity of 0.98, because the majority of 
the stripped mass passes within 100 pc of the galaxy centre 
and is therefore removed from the simulation. By compar- 
ing these profiles with those in Figure 3, it can be seen that 
the result of allowing the gas particles to interact with each 
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other in a more realistic manner is to broaden the distribu- 
tions, and to move the innermost gas particles inwards by 
about a factor of two in radius. Some of this difference is 
due to the more realistic treatment of the stripping process, 
some is due to treating the gas dynamics more realistically 
and some is due to the artifical viscous effects introduced 
by our approximate scheme for treating the gas dynamics. 
We believe that the viscous effects are not significant here. 
The mean radius of the final gas distribution is within 20% 
of the results obtained using the simple prescription of the 
previous section, which implies the total energy loss is simi- 
lar. This approximate similarity between the two approaches 
serves to indicate that any viscous effects introduced due by 
our grid-based dissipation scheme are reasonably small. 

3.3 High energy orbits 

The above results make it clear that in order to form a dust 
disc of the size of a few hundred parsecs of the kind which is 
seen in NGC 3379, the orbit of the hypothesized incoming 
satellite galaxy is not of the form of a low energy, low eccen- 
tricity orbit which goes through multiple pericentre passages 
before the satellite is disrupted. Rather the initial orbit of 
the incoming satellite needs to have a sufficiently high ec- 
centricity (i.e. low angular momentum) that the satellite is 
likely to be completely disrupted on its first pass. 

For a more detailed investigation of the high eccentric- 
ity regime, it is more useful to parameterize the problem in 
terms of an impact parameter b and an initial velocity Vb- 
In addition, to model the approach of the satellite galaxy 
to the host's core, we now start the satellite's orbit at a 
galactocentric radius of 80 kpc. We use the same procedure 
as outlined above (Section 3.2) to treat the gas stripping 
and subsequent gas dynamics, but, since the pericentric dis- 
tances are now considerably reduced, we now decrease the 
inner radius at which particles are removed from the com- 
putation from 100 to 10 pc. This is necessary because of the 
change in the nature of the orbits of the gas particles. When 
we were considering low energy orbits, the particles which 
reached 100 pc were on essentially circular orbits and are 
assumed to become part of the nuclear gas disc, but in the 
high eccentricity case, particles reach 100 pc with velocities 
much greater than circular and as such do not necessarily 
get incorporated into the disc. 

With these initial conditions, we find that a disc of sim- 
ilar size to that observed in NGC 3379 can be reproduced 
if the specific angular momentum of the incoming satellite, 
bVo is around 50 kpc km s _1 . The surface density profile of 
the disc obtained using this value for the orbital angular 
momentum is shown in Figure 5. In this situation, around 
1% of the gas particles fall within the radius 10 pc and so 
are removed from the calculation. 

We conclude that in order to obtain a dust disc the ap- 
proximate size of the one observed it is necessary for the 
satellite galaxy to be on an orbit which is almost precisely 
oriented towards the galactic nucleus. In order to estimate 
of the likelihood of such an occurrence, it is of course nec- 
essary to have some knowledge of the space and velocity 
distribution of the prospective mergee population. This we 
do not have. Numerical computations of galaxy formation 
through mergers should in principle be able to give some in- 
sight into this (e.g. Steinmetz 1999), but current simulations 



do not currently have high enough mass resolution for our 
purposes. 

We may, however, use the above results to estimate 
the fraction of those merging satellites disrupted on their 
first pericentre passage which give rise to a dust disc of the 
kind observed in NGC 3379. If we assume that the satellite 
galaxy, as it falls into the main galaxy, has a velocity ap- 
proximately equal to that of the dispersion velocity in the 
main galaxy, which is about 200 km s~\ then for this satel- 
lite give rise to the observed gas disc, we have seen that we 
require bVo ~ 50 kpc km s _1 , and thus that the pericentre 
of the orbit to be b ~ 250 pc. For the satellite to be dis- 
rupted on the first pass through the main galaxy, we have 
seen (Section 2.3) that we require the pericentre distance to 
be less than around b ~ 1 kpc. If we assume that incoming 
satellites have their impact parameters evenly distributed 
over this range then this would imply that only of the or- 
der of 4% of those satellites disrupted on the first pass are 
on appropriate orbits to form a disc such as observed in 
NGC 3379. Of course, the majority of satellites have larger 
impact parameters. They would be captured by the main 
galaxy, and would then slowly work their way inwards in 
the manner described in Section 2.3, distributing their gas 
and dust widely through the host galaxy. This should lead 
to observable effects (Section 5). This in any case we ex- 
pect that for every merging satellite which could produce 
the observed disc, there are many more which could not. 



4 MERGING WITH A SPIRAL GALAXY 

For the reasons discussed in the Introduction, we now extend 
the foregoing analysis to consider the merging of a small 
satellite galaxy with a spiral host galaxy. Our interest in this 
case centres not so much on where the stripped gas may end 
up (since it likely in any case to collide with and merge with 
the gas already in the gas disc of the spiral), but rather on 
the fate of the nucleus of the satellite, or more especially, on 
the fate of the black hole presumably contained within it. 



4.1 Galactic models 

4-1.1 The spiral 

We choose to use a simple model of a typical disc galaxy, 
in which the galaxy is modelled using four axisymmetric 
components (representing the central black hole/star clus- 
ter, the bulge, the disc and the dark halo), each of which 
is represented by a Miyamoto-Nagai potential (Binney and 
Tremaine 1987, Equation 2-50a), which takes the form 

$(R, *) = ; GM (13) 



sJr 2 + (o + Vz 2 + fc 2 ) 2 



where we use a cylindrical coordinate system (R,(f>,z). 
The density due to this potential is 

,„ , b 2 MaR 2 + (a + 3^/PTT 2 ") (a + V6H?) 2 , 

p(R>*) = 2 , 5/2 M14) 

(R 2 + (a + Vb 2 +z 2 ) J (b 2 + z 2 f 2 

We note that this reduces to a Plummer sphere in the 
case of a = 0. From these expressions, and using the as- 
sumption that the internal galactic dynamics consists of an 
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Figure 4. The surface density distribution of gas stripped according to the model in Section 3 for orbits with initial eccentricities of 
0.5 (left panel) and 0.75 (right panel). Results for an eccentricity of 0.98 are not shown as the majority of the mass is removed from the 
simulation; see Section 3.2 for details. 
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Figure 5. A plot of surface density against radius for the disc formed when bVo is 50kpc km s 1 (left panel) and 100 kpc km s 1 (right 
panel). 



isotropic velocity dispersion coupled with a mean azimuthal 
streaming motion, we can calculate the velocity dispersion 
and mean streaming velocity of each component, via (Bin- 
ney and Tremaine, 1987, Equations 4-65 and 4-66): 



a ( ,z) = ~ J Kz) ] p( z) — dP — ( } 



vl(R,z) = R 



+ 



p(R, 
d$(R,z) 



8R 
R d 



p(R,z) OR 



, d$(R,z') , 
P(R,z) dz , dz (16) 



As before, enlarging the parameter space to make more 
complicated assumptions about the stellar distribution func- 
tion is not warranted given the other simplifications we 
adopt here. In any case, the results are not likely to be 
strongly affected unless the stellar distribution function is 
extremely anisotropic. 

For our disc galaxy, we take the parameters used by So- 
me (1996) for the Milky Way which assumes: (i) a spherically 
symmetric (a = 0) nuclear star cluster of mass 5 x 10 9 M , 
and scale b = 120 pc (ii) a spherically symmetric bulge of 
mass 10 10 M and scale b — 750 pc; (iii) an axisymmetric 
disc of mass 1.6 x 10 11 Mq , radial scale a = 6 kpc and verti- 



cal scale b = 500 pc; and (iv) a spherically symmetric dark 
halo of mass 3 x 10 11 M Q and scale b — 15 kpc. 



4.1.2 The satellite 

As above, we model the satellite galaxy in this case as a 
Plummer sphere, subject to tidal stripping but with the ad- 
dition of a central point mass to represent a nuclear black 
hole and surrounding star cluster. We increase the mass of 
the satellite to be 4 x 10 9 Mq but retain the same scale 
length of b — 0.4 kpc. This is consistent with the models 
of Sagittarius of Helmi and White (2001). This represents 
the inclusion of dark matter in the model for the dwarf. 
This could have be added into the previous simulations, but 
would have had little effect as the satellite was very rapidly 
destroyed by tidal stripping. We set the central point mass 
to be 1 x 10 7 Mq , although we shall see later that our sim- 
ulations are basically independent of this. 



4.2 Orbital dynamics 

We model the orbital dynamics of the satellite in a manner 
using basically the same ideas as used above (Section 2.3). 
The satellite galaxy is subject to the gravitational influence 
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of the main galaxy, to dynamical friction and to tidal strip- 
ping. 

The gravitational force due to the main galaxy is 
straightforward to calculate. Due to the axisymmetry of our 
model, there are only two independent components to con- 
sider. The radial force component (in a cylindrical sense, 
that parallel to the plane of the disc) is 



Fr 



9$ 
dR 



GMR 



3/2 



and the vertical component (perpendicular to the disc), 



(17) 



F - 9<E> - 



GMz (a + Vb 2 + z 2 ) 



(i? 2 + (a + Vb' 2 + z 2 ) 2 ^j ' i/2 (a + Vb 2 + z 2 ) 



.(18) 



The force in both directions is calculated for each of 
the four components of our model galaxy and these forces 
summed to produce the resultant force on the satellite due 
to the main galaxy. 

To calculate the dynamical friction on the satellite, we 
first calculate the mean streaming velocity, v^, of the disc 
component of the galaxy at the satellite's current location 
using Equation 16 (note that all other components have zero 
mean streaming velocity since they are spherically symmet- 
ric) and from this we calculate the relative velocity of the 
satellite which respect to the disc. Second we compute the 
local velocity dispersion of each component of the galaxy. 
Then, finally, we use the standard gravitational drag for- 
mula (Equation 7) to calculate the dynamical friction force 
on the satellite. Once again, the contribution from all four 
components of the galaxy is then summed, and this added 
to the gravitational force to produce the total force on the 
satellite. 

After each timestep in our integration, we apply simple 
tidal stripping to the satellite. We calculate the mass internal 
to the satellite's current position, via 



(19) 



f"\J R 2 +z 2 pt\ 

M(R,z) = 2ir / p{r,z)r 2 sin 6&6Ar 

Jo Jo 

and then use this mass to calculate the current tidal radius 
of the satellite, as given in Equation 10. Whilst this is not 
strictly accurate, it is an adequate approximation for our 
current purposes. As we shall see below, the stripping of the 
satellite is not the most important factor in these simula- 
tions. 

Again, we use a variant of the 'odeint' routine from 
Numerical Recipes (1992) to control the timestep in the sim- 
ulation, so that we use small timesteps when necessary, but 
larger ones when possible to speed up our simulation. 



4.3 Results 

Since we know from our previous set of calculations that the 
only satellite trajectories which are likely to have a chance 
of reaching the centre of the host galaxy are those which 
are almost radial, we assume that the satellite galaxy is 
initially approaching the host galaxy from a large distance 
(> 100 kpc). 

To be specific, we define the centre of the Seyfert as 
the origin of our coordinate system, and then consider the 



trajectories of small, incoming galaxies which start on the 
plane which has a normal vector at an angle 90 — a from 
the axis of rotation of the Seyfert and has a point of clos- 
est approach of 100 kpc to the origin. We consider satellite 
orbits which start in this plane, and which have their initial 
velocity being perpendicular to the plane. We shall call the 
angle a the 'approach angle', as it corresponds to the angle 
between the initial velocity vector and the plane of the host 
galaxy. Thus a — 90° corresponds to an orbit which is ini- 
tially directed at right angles to the disc of the host galaxy. 
For this paper, we investigated orbits with approach angles 
of 0° , 15° , 30° , 45° , 60° , 75° , 85° , and 90° and with two sets 
of initial velocities: 100 kms -1 and 200 kms -1 . 

We continue to track the satellite's orbit either until the 
satellite has settled into the plane of the disc, after which 
orbital evolution is very slow, or until the satellite reaches 
a radius, r sp h at which the potential of the host galaxy is 
spherical and so there will be no further changes in the direc- 
tion of satellite's orbital angular momentum. For our model, 
we take r sp h to be 200 pc. 

The main outcome of our orbital computations is that 
most satellites settle fairly quickly into the orbital plane 
of the disc galaxy. This comes about simply because the 
gravitational drag due to the plane is relatively large because 
the dispersion velocity of the disc component is relatively 
small. The fraction of the initial conditions which result in 
the satellite reaching the central regions of the main galaxy 
before it settles into the plane of the disc is relatively small, 
but not vanishingly so. Those that do reach the centre are 
generally those for which their first dynamically significant 
impact with the disc is in the centre of the galaxy. Three 
orbits, chosen such that two of them (the solid and dotted 
lines) are those which reach the centre of the galaxy, are 
shown in Figure 6. Once the satellite has settled into the 
plane its orbital evolution becomes very slow, and does not 
in general reach the centre within a Hubble time. We should 
note that for our approximations of a static potential for 
the disc to be valid, we require that the satellite not be 
more massive than the mass it interacts with when passing 
through the disc. The satellite will interact with material 
in a cylinder of radius of order of GM/v 2 , where v is the 
velocity of the satellite when it passes through the disc. This 
mass is always greater than the current mass of the satellite 
when passing through the disc, except for the earliest passes 
through the disc when there is no significant drag. 

For each orbit which reached r sp h before settling into 
the plane of the disc, (i.e. which reached the centre in less 
than a Hubble time) we recorded the angle, /3, which the 
satellite's final orbital angular momentum vector made with 
the normal to the plane of the disc. Thus /3 = 0° means 
the satellite's angular momentum is parallel to that of the 
rotational angular momentum of the disc, i.e the final orbit 
of the satellite lies in the plane of the galaxy disc. 

In general, the plane (a — const.), from which we 
launch our initial trajectories contains a small number of 
regions (typically around 5) from which the satellite reaches 
the centre of the host galaxy, but these are spread out over 
a relatively large area. To find the cross-sectional area A (a) 
at each approach angle, we first undertook a low-resolution 
scan (out to a radius where satellite simply flew past the 
galaxy), to find the regions in which the satellite reached 
the centre before reaching the plane of the galaxy. We picked 
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Figure 6. Three orbits in the disc galaxy; the disc is perpen- 
dicular to the page. The approach angle, a is 90° and the initial 
velocity is 100 km s -1 in all three cases. 



Figure 7. The contribution to the cross-sectional area for the 
satellite reaching r sp h before settling into the disc of the host 
galaxy for various approach angles a (see Section 4.3). 



out the areas of interest by eye, and then investigated those 
regions at higher resolution. The cross-sectional area was 
then calculated simply by summing the area of these re- 
gions. The exception to this is the face-on case, which is ro- 
tationally symmetric in the plane of the disc. Here we simply 
undertook a high resolution scan along one line parallel to 
the plane of the disc and then calculated the cross-sectional 
area from the circular annuli produced. In Figure 7, we plot 
the contribution A (a) cos a to the total cross-sectional area 
made by each approach angle a. The factor of cos a accounts 
for the fact that we expect fewer satellites to have approach 
angles close to 90° , assuming the incoming flux of satellites 
is isotropic. 

We have also determined the angle (3 between the or- 
bital angular momentum of the satellite's orbit when it 
reaches the centre (r < r sp h) and the normal to the plane of 
the galaxy disc. Figure 8 shows the distribution of /3 for dif- 
ferent approach angles a. We note that there is little change 
with a. The absence of orbits which end with a low value of 
P is due to the fact that these orbits are close to the plane 
of the disc, and so quickly settle into the plane of the disc 
and are never reach the radius r sp h. 

We also compute an overall distribution for the angle f3 
for those mergees which reach the centre of the host galaxy. 
In doing so we encounter the problem that the distribution 
of the initial orbits of the satellites' orbits is unknown. In the 
absence of any strong evidence to the contrary, we assume 
that the incoming merging satellites approach from random 
directions. With this assumption, we have computed the dis- 
tribution for [3, and this is as shown in Figure 9. We note 



that this distribution is reasonably uniform in /3, except for 
a zone of avoidance for angles /3 < 20° . 

We have repeated the above experiments with the black 
hole in the centre of the dwarf galaxy having a mass from 
10 5 M to 10 8 M (as opposed to the 10 7 M used previ- 
ously), and obtain no significant differences. This may have 
been expected from the fact that the satellite is not stripped 
down to a small mass by the time it reaches a radius of 200 pc 
from the centre of the Seyfert. 

4.4 Comparision with observations 

A comparision of our model with the available observational 
data is not entirely trivial. In order to obtain a quantitative 
result, we follow the method of Kinney et al (2000, and ref- 
erences therein). A brief description of the method is given 
here. 

Two parameters are observed for every galaxy used in 
Kinney et al.: i, the inclination of the galaxy's disc to the 
plane of the sky, and S, the difference between the position 
angles of the galaxy's major axis and that of the jet when 
projected onto the plane of the sky. Each (z, 8) pair con- 
strains the jet to lie on a pair of great circles on a sphere 
about the nucleus of the galaxy. The angle /3 between the jet 
and the rotation axis of the disc varies around these great 
circles. We assume that the chance of a jet being at any 
point on either of these great circles is given by the proba- 
bility of our simulations producing the value of /3 (as shown 
in Figure 9) for that point on the great circle. At a fixed in- 
clination, different values of 8 give great circles which have 
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Figure 8. The distribution of the angle /3 between the angular 
momentum vectors of the infalling satellite and the central black 
hole of the Seyfert galaxy for different approach angles a of the 
satellite's initial orbit: 0°(dashed), 30°(dotted) and 60°(solid). 
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Figure 9. The overall distribution of the angle /3 assuming that 
the initial distribution of orbital inclinations is uniform. 



different orientations in space, and hence we calculate the 
probability of obtaining a given value of delta given the in- 
clination, p(8\i). We then integrate this distribution from 
8 = to 8 = 6k where <5fc is the observed value of 8 for the 
galaxy k to obtain the centile point Ck- The distribution of 
centile points across the data set should then be uniform 
over the range [0,1]. Full details of this technique can be 
found in Kinney et al. (2000) . 

We perform a K-S test to compare our model with the 
the 60 fim data of Kinney et al. (2000) . We find that the 
data is consistent with being drawn from our model at the 
97% level. 



5 DISCUSSION 

We have undertaken an investigation into the possibility 
that the activity of, and jet direction in, at least some active 
galactic nuclei might be fuelled by the infall and merger of 
small satellite galaxies. 

5.1 Cross-sections for mergers 

We considered first the evidence that many of the host galax- 
ies active nuclei in the 3CR catalogue are found to have discs 
of dust around their nuclei. Taking this as prima facie evi- 
dence that the nuclear activity might be related to infall, we 
have modelled the disruption of a small satellite galaxy by 
a larger elliptical, taking the particular case of NGC 3379, 
which while not in the 3CR catalogue, does show a radio 
bright nucleus and does display a dust disc. 

Our modelling procedure is of necessity somewhat 
crude. We have adopted a realistic model for the inner re- 
gions of the host elliptical where the disruption of the in- 
falling satellite takes place, but have not taken proper ac- 
count of a possible massive, but low density, dark halo. We 
considered idealised models for the internal structure of the 
satellite, one with and one without a central density cusp, 
but found that the results did not depend significantly upon 
which models we used. Gravitational drag on the orbit of 
the satellite was modelled approximately using the standard 
Chandrasekhar formula, suitably parameterized by compar- 
ison with full numerical computations taken from the liter- 
ature. 

Probably the least realistic part of our modelling proce- 
dure is the description of the disruption of the satellite and 
the subsequent dynamics of the gas it contained. For exam- 
ple, we have assumed that the main galaxy is sufficiently 
devoid of gas that ram pressure stripping of gas from the 
satellite can be ignored, and that the effect of the interac- 
tions between the stripped gas and the gas already present 
may be ignored. This may not always be the case. Even so, 
because our description does contain much of the relevant 
physics it is unlikely that our results will be grossly in error. 
The infalling satellite is disrupted at the point where tidal 
effects of the host come into play using simple truncation of 
the underlying model. This ignores the dynamical response 
of the satellite to the truncation process, and so, if anything, 
tends to overestimate the radius at which disruption occurs. 
On the other hand, we have implicitly assumed that the gas 
and stars in the satellite occupy the same volume, and thus 
strip gas particles in proportion to the mass removed. If the 
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gas were more centrally located than the bulk of the mass, 
then our procedure would underestimate the radius at which 
the gas from the satellite is deposited. However, as can be 
seen from Figure 2, the radius at which tidal disruption takes 
place is fairly well-defined, with the radius by which 90 per 
cent of the satellite has been tidally removed being almost 
independent of the model for the internal structure. The 
only way to bring about a significant change in the radius 
at which the gas is deposited is to use the cusped model and 
to assume that all the gas lies within the radius occupied 
by the central one percent of the total galaxy mass. This 
would imply for a gas to mass ratio of a percent that the 
gas was dynamically important in the nuclear region of the 
satellite. The modelling of the dynamics of the gas once it 
has been stripped conserves mass and angular momentum, 
but uses a simple algorithm to allow intersecting streams of 
gas to shock and lose energy by radiation. This is probably 
an adequate approximation for the current paper. 

Despite the relative crudeness of our modelling proce- 
dure, we are able to draw some conclusions which we believe 
are relatively robust. We find that while it is possible to re- 
produce the dust disc structure seen in NGC 3379 via the 
infall of a satellite galaxy, it seems likely that this type of 
event is very rare. Our initial simple model for the satel- 
lite, which started within the galactic halo on an orbit with 
a given eccentricity, was not able to reproduce the struc- 
ture seen in NGC 3379. A more realistic model, in which 
the satellite starts well outside the galaxy with a small im- 
pact parameter was found to be capable of reproducing the 
observed structure, provided that the satellite galaxy has 
an initial angular momentum about the galactic nucleus of 
around 50kpc km s _1 . Given the lack of observational or 
theoretical evidence as to the distribution of the orbits of 
dwarf satellites, we cannot make a quantitative statement 
about how likely this is, but we do not believe that this 
event is common. If we assume a typical velocity of around 
200 km s -1 for the satellite when at infinity, this means the 
initial orbit must have an impact parameter of the order 
of 250 pc, which would represent a small fraction of phase 
space if the initial orbits were uniformly distributed. Even if 
the orbital distribution were moderately biased towards low 
angular momentum orbits, this would still imply that the 
formation of a structure similar to that in NGC 3379 being 
formed by this method is a rare event. 

With regard to powering a single AGN, this model pro- 
duces results which are not inconsistent with current theo- 
ries. To order of magnitude level, if an event forming a disc 
similar to that observed in NGC 3379 occurs, we find that 
a few percent of the gas particles reach a distance of 10 pc 
from the nucleus. At this point it is reasonable to assume 
that they will be swallowed by the central black hole within 
a reasonably short time span (Shlosman et al. 1990; Bekki 
2000). If we assume that a few percent of the mass of the 
original dwarf galaxy was gas, this gives around 10 6 Mq be- 
ing swallowed by the black hole, which, at an accretion rate 
of 0.1 M Q /yr, will power a low level AGN for around 10 7 yr. 
This is not unreasonable. 

However, if we further assume that all galaxy activity is 
triggered only by events of the above type, then the model 
becomes harder to sustain. If we assume that every merger 
triggers nuclear activity which lasts around 10 7 years, and 
if around 1% of elliptical galaxies are active, then this im- 



plies that each elliptical galaxy has undergone about 10 such 
events over the lifetime of the universe. However, if only 
a few percent of mergers occur at sufficiently low angular 
momentum to give rise to nuclear activity, then the total 
number of satellites swallowed by the host over its lifetime 
rises to around 300 to 1000. If each satellite has a mass of 
around 10 9 M as we assume here, this implies the addi- 
tion of a substantial amount of mass. Even if our estimates 
are too pessimistic with regard to the efficiency with which 
merger process is able to fuel the nucleus, it is clear that 
we would nevertheless expect that essentially all elliptical 
galaxies should show some evidence of such a merger hav- 
ing occurred, and perhaps more than one. In this regard, we 
note that whilst van Dokkum and Franx (1995) found dust 
structures in approximately half of their sample, more recent 
more sensitive surveys (Pastoriza et al. 2000) are finding 
dust in a higher proportion of elliptical galaxies, including 
some in which dust remained undetected in the van Dokkum 
and Franx survey. 

5.2 Directionality 

Having established that the likelihood of any particular 
merging satellite galaxy reaching the host nucleus is rela- 
tively small, we then turn our attention to incoming orbits 
for which the successful mergees actually reach the nucleus. 
The motivation for this part of the investigation is the dis- 
covery (Kinney et al. 2000) that the jet directions in the 
centres of Seyfert galaxies appear to bear no relation to the 
orientation of the discs of the host galaxies. The possibility 
we are investigating here is whether the mergees reaching the 
centre do so with an angular momentum distribution which 
could produce the observational data observed by Kinney et 
al. 

For this case we have taken the host galaxy to be a spi- 
ral galaxy, and have employed the same procedures for or- 
bital decay and tidal stripping as before. Again we find that 
the cross-sectional areas for merging of the satellite with 
the host nucleus are relatively small, but we now concen- 
trate our analysis on the angular momentum of the orbits of 
those satellites (or what is left of them) when they reach the 
nuclear regions of the host. Those satellite nuclei which do 
reach the centre in a finite time (that is, before settling into 
the disc) have the angle between their angular momentum 
vector and the disc normal (/3) distributed approximately 
randomly, as is shown in Figure 9. The deficit at low values 
of p comes about because these satellites are close to, or 
in, the disc of the Seyfert galaxy, and as such rapidly settle 
down into the plane. We have noted that is distribution is 
in good agreement with the data presented in Kinney et al. 
(2000). 

The simulations we have performed here give us a dis- 
tribution for the angle /3. However, in order to relate these 
results to the observed jet distribution, we need the resul- 
tant spin of the black hole formed between the merger of 
the two holes. For the general case, this is an unsolved prob- 
lem in general relativity and hence we again adopt a simple 
prescription. We assume the mass of the dwarf's black hole, 
Mn, is much less than the mass of the Seyfert's black hole, 
Magn and can then use the prescription from Colbert and 
Wilson (1995) to first order to obtain the spin of the daugh- 
ter of the merger as 
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L S = 2 ^ GMoMAGN (20) 
c 

where L$ is the spin of the resultant hole. This applies 
only for a merger between two black holes which were not 
originally spinning, which we note is not the case here. If 
Md/Magn = 0.1, then the resultant hole is a Kerr black 
hole with spin parameter a — 0.35, where a is the angular 
momentum of the hole in units of GM 2 /c and lies between 
and 1. Thus the accumulation of a few satellite nuclear 
black holes by the black hole in the AGN nucleus, from or- 
bital directions which are more of less randomly distributed 
can lead to random orientations for the spins of the AGN 
nuclear black holes. A more detailed discussion of this is 
given by Merritt (2002). 



6 CONCLUSIONS 

We have investigated the dynamics of the merging process 
in the minor merger hypothesis for active galactic nuclei. 
Our analysis has resulted in two main conclusions. 

First, for the satellite galaxy to be able to merge directly 
with the nucleus of the host galaxy (for example, to give rise 
to the compact dust discs which are seen in early type active 
galaxies) requires the initial orbit of the satellite to be well 
aimed. The corollary of this is that for each merging satellite 
which gives rise to a compact nuclear dust disc, there must 
be many which do not. Similarly, in disc galaxies, for every 
minor merger which causes an active nucleus, there must 
be many mergers which simply give rise to a merger of the 
incoming satellite with the galaxy disc. 

Second, in the case of the host galaxy being a disc 
galaxy, if the initial orbits of the satellites are randomly 
oriented with respect to the host galaxy, then the orbits of 
those which reach the host nuclear regions in a reasonable 
time, are also fairly randomly oriented once they reach the 
nucleus. We note that this result might be able to provide an 
explanation of why the jet directions in the nuclei of Seyfert 
galaxies are apparently unrelated to the plane of the galaxy 
disc. The distribution of angles (3 between the final orbit of 
the successful merger candidates and the Seyfert disc which 
we find once we have averaged over an assumed isotropic 
incoming satellite distribution is consistent with the obser- 
vational data presented by Kinney et al (2000). 
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